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We consider the initial-boundary value problem for systems of quasilinear wave equations on do- 
mains of the form [0, T] x E, where S is a compact manifold with smooth boundaries 9E. By using 
an appropriate reduction to a first order symmetric hyperbolic system with maximal dissipative 
Qy ' boundary conditions, well posedness of such problems is established for a large class of boundary 

conditions on 9S. We show that our class of boundary conditions is sufficiently general to allow for 
a well posed formulation for different wave problems in the presence of constraints and artificial, 
nonreflecting boundaries, including Maxwell's equations in the Lorentz gauge and Einstein's gravi- 
tational equations in harmonic coordinates. Our results should also be useful for obtaining stable 
finite-difference discretizations for such problems. 
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INTRODUCTION AND MAIN RESULTS 



Motivated in part by the numerical computation of spacetimes on a finite domain with artificial boundaries, the 

initial-boundary value problem (IBVP) in general relativity has started to receive a lot of attention during the last few 

^_i years (see [l| for a review). A well posed IBVP for Einstein's vacuum field equations was formulated for the first time 

__j' by Friedrich and Nagy [3] based on tetrad fields and the theory of quasilinear, symmetric hyperbolic systems with 

maximal dissipative boundary conditions [3|, y, |5| . More recently, Kreiss and Winicour [6] formulated a well posed 

^""i ' IBVP for the harmonic gauge formulation of the Einstein vacuum equations which casts the field equations into a set 

^^ , of ten coupled quasilinear wave equations subject to four constraints. There are two key ideas behind the result of 

L«. ' [6|. The first one is the realization that the wave equations, when viewed as first order pseudodifferential equations, 

-vj ] has a non-characteristic boundary matrix. This allows application of the boundary value theory for such systems 

r-o ■ developed by Kreiss in the 1970's Q- Second, the formulation of boundary conditions for the frozen coefficient form 

• ■ of the harmonic Einstein equations which ensure constraint propagation and satisfy the estimates required by the 

\-^ ,: Kreiss theory. The well posedness of the system and the generalization to the quasilinear case can then be established 

QQ ■ using the theory of pseudodifferential operators (see, for instance, ^]). 

f~^ , In a subsequent paper [9|, similar results were obtained via more mundane energy estimates which follow by 

L* ' integration by parts, without resort to the pseudodifferential calculus. For this, a non-standard energy norm is 

. !^ , constructed which is based upon the choice of a particular time-like direction adapted to the boundary conditions 

Si^ ■ being imposed. With respect to this energy the Kreiss- Winicour boundary conditions are maximally dissipative and 

?H ' so standard well posedness theorems apply even in the quasilinear case [5|, llOj . Besides being a simpler proof, or 

at least a proof that can be followed completely by a reader not familiar with the pseudodifferential techniques, it 

implies similar results for the stability of finite difference approximations to Einstein's equations in the harmonic 

gauge. This follows from considering the semidiscrete system of ordinary differential equations in time obtained by 

substituting finite differences for spatial derivatives. If the semidiscrete system is stable, then for appropriate time 

discretizations the fully discrete system is guaranteed to be stable [ll| . The stability of the semidiscrete system can be 

established by the use of finite difference operators satisfying summation by parts [I4I , the counterpart of integration 

by parts, by mimicking the steps leading to the continuum energy estimate. A summation by parts algorithm for 

the harmonic Einstein IBVP was developed for homogeneous boundary conditions 13^ and verified to be stable in 

numerical tests [14] . The results of [y| provide a means to prove stability in the inhomogeneous case. 

In this paper we present a more general and geometric version of the foregoing results which applies to coupled 
systems of quasilinear wave equations with a certain class of boundary conditions. The well posedness of the resulting 
IBVP is established by reducing the wave system to first order symmetric hyperbolic equations subject to maximally 
dissipative boundary conditions. As we show, our class of boundary conditions is sufficiently fiexible for obtaining well 
posed IBVP formulations for different models of isolated systems in physics, including the wave equation. Maxwell's 
equations and the Einstein field equations. 
In what follows we present the main results. 



A. Main theorem 

Let T > 0, and denote by S a d-dimensional compact manifold with smooth boundaries d'E. The type of system 
our results apply to is a set of quasilinear wave equations on M = [0, T] x E coupled both by lower order terms and 
in the principal part, by a change in the characteristic directions via a metric which can depend on the local value 
of the fields involved. More precisely, let tt : i? ^ M be a vector bundle over M with fibre R^, let Va be a fixed, 
given connection on E and let gat — Qabi^) be a Lorentz metric on M with inverse g°'^{^) which depends pointwise 
and smoothly on a set of fields $ = {^^}a=i,2,...n parameterizing a local section of E. Our signature convention 
for gab is (— ,+, •••,+)■ We shall also assume that each time-slice E* — {t} x E is space-like and that the boundary 
T = [0, T] X 9E is time-like with respect to gab{^)- In the following, we will refer to local sections in E as vector- valued 
functions over M. We will also assume the existence of a positive-definite fibre metric Hab on E. We consider a 
system of quasilinear wave equations of the form 

5''''($)VaV6$^ = 5^($,V$), (f) 

where 5 ($,V$) is a vector- valued function which depends pointwise and smoothly on its arguments. The wave 
system ([T]) is subject to the initial conditions 

where $g and Ilg are given vector- valued functions on Eq, and where n^ — n''($) denotes the future-directed unit 
normal to Eg with respect to gab- In order to describe the boundary conditions, let T" = T'^{p, <&) be a future-directed 
vector field on T which is normalized with respect to gab and let N'^ = N"-{p, $) be the unit outward normal to T 
with respect to the metric gab- We consider boundary conditions on T of the following form^: 

[T^ + aN''] wM\r = c'^^b Va$^|^ + d^B <^''\r + G^, (3) 

where a — a(p, $) > is a strictly positive, smooth function, G — G (p) is a given, vector- valued function on 
T and the matrix coefficients d^^B = c"'^b(p, $) and d-^B — d^B{p,^) are smooth functions of their arguments. 
Furthermore, we assume that c"' b can be made arbitrarily small in the following sense: Given a local trivialization 
(fi : U X M.^ 1-^ TT^^{U) of E such that U C M is compact and contains a portion U of the boundary T, and given 
£ > 0, there exists a smooth map J : U —t GL{N,'R),p i-^ {J^b{p)) such that the transformed matrix coefficients 

C B -— J CC D [J ) B 

satisfy the condition 

hAB£'''c{<^)c'"'Di'^)Va^Vb'' < ehABe''\<i>)Va^Vb'' , (4) 

for all vector- valued one- forms Va on U, where here and in the following, Cab refers to the Euclidean metric Cab = 
gab + '^TaTf, which is defined for points on T. 
The main result of this paper is: 



Theorem 1 The IBVP U\2\3\) is well posed. Given T > and sufficiently small and smooth initial and boundary 
data $0 , IIq and G satisfying the usual compatibility conditions at 9Eo, there exists a unique smooth solution on 
M satisfying the evolution equation {Ip, the initial condition (0j and the boundary condition (0). Furthermore, the 
solution depends continuously on the initial and boundary data. 

A common situation in which the condition ^ is automatically satisfied is given in the following 

Lemma 1 Let U d T be an open subset of T such that U is compact. Assume there exists a smooth map J : U —i- 

GL{N,M.),p i—> (J sip)) over U such that the transformed matrix coefficients c"" b '■= J cc"'^D {J~^) B are in 
upper triangular form with zeroes on the diagonal, that is 



d'''^B = 0. B<A. 



Then, the condition Ol) is satisfied on hi. 



^ We adopt the Einstein summation convention for the lower case Latin abstract spacetime indices a, b, c, ... as well as for the Capital 
indices A, B, C, ... on the fibre of E. 



Proof, (cf. The proof of the Liapunov stability theorem) In order to simphfy the notation we use a matrix notation 
and write c" = Jc°-J-\ Let 5 > 0, and define Dg := diag{l, 6,6^, ■■■,S^-^) and Js := DJ^J. Then, c^ := Jgc^J^^ = 
DJ^c°'Ds has the components {c'g)"^B = S^^'^c°'^b, where here, 6^~"^ refers to the {B — A)'th power of S. Since 
£"■^3 = for _B < A we have Cg — 0{S), and Cg satisfies the condition ((4|) provided (5 > is chosen small enough, q 

The proof of theorem [1] is given in sections |TI] and IIIII In order to illustrate the ideas on a simpler example, we 
start in Sect. ITTlwith the wave equation on a fixed background metric gab, and analyze the general case in Sect. IIIII 

Since many physical systems can be described by systems of wave equations, theorem [T] should have many applica- 
tions. In the following, we mention two such applications for the initial-boundary value formulation of isolated systems 
with constraints. The physical motivation for the choice of nonreflecting boundary conditions in these examples is 
described in detail in section Hvl 

B. Maxwell's equations in the Lorentz gauge 

The first application describes an electromagnetic field on the manifold M — [0, T] x S with a fixed background 
metric gab and corresponding Levi-Civita connection Vq. As before, we assume that each time-slice Et ~ {t} x S 
is space-like and that the boundary T = [0, T] x 9S is time-like. In the Lorentz gauge C := VbA'' — 0, where A'' 
denotes the 4- vector potential. Maxwell's equations assume the form of a system of wave equations, 

g^'S/a^bA- = R'dA' - r , (5) 

where Rab denotes the Ricci tensor belonging to the metric gab and J"^ is the four-current. ([3]) implies that the 
constraint variable C obeys the following equation 

g^'WaVbC = - VV, . (6) 

Therefore, the imposition of the boundary condition C\j- — and the satisfaction of the continuity equation VJc = 
imply that any smooth enough solution of ([5]) with initial data satisfying 

C|so = 0, n^VaC\^, = 0, 

satisfies the constraint C — on M since in this case the constraint propagation system ([6]) is homogeneous. 

Asymptotically nonreflecting boundary conditions at T = [0, T] x E, in the sense of II VI can be formulated by first 
introducing a null tetrad {K°',L°',Q°',Q'^} which is adapted to the boundary. Let T" be a future-directed time-like 
vector field tangent to T normalized such that gabT'^T'' — —1, let TV" denote the unit outward normal to T with 
respect to gab and complete T" and N"' to an orthonormal basis {T°, A^°, V"", W^} of TpM at each point p E T. Then, 
we define the null vectors 

K" := T" + N" , L" -.^ T" ~ N" , Q" -.^ V" + iW, Q" := V - iW, 

where i — \/— T. These vectors may be smoothly continued in a small region inside the domain, for example by parallel 
transport along the normal direction to the boundary. In this way, one obtains a local null basis of TM. Finally, let 
r denote the areal radius of the cross sections d'Et . This function can also be continued in a small region inside the 
domain by parallel transporting dT,t along the normal direction. The following boundary conditions are derived in 
section HYB] 

^K-KbX/air^A") 
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{K^Qb - Q''Kb)VaA\ = qg , (8) 

{K''Lb + L''KbVaA''-Q''Qb-QaQ'')^aA\ = 0, (9) 

where qk and qg are given real and complex scalars on T. The first condition is a gauge condition, the second condition 
controls the electromagnetic radiation through T and the third condition enforces the constraint C = g^^VaAf, ~ 
onT. 

The evolution equation ([5]) has the form ([1]) where E is the tangent bundle over M, and the boundary conditions 
d7l8l9)) have the form © with 



a — 1, 

G^ = ^ [-L^qx + Q'qg + Q'qq] 



2g(°Q^)i^rf + L^K'Kd - K' {Q^Qd + Q^Qd) , d'd = (A'^V^ log r)i^i^rf 



Since 

ffc^j^d ^ _iaj^c _ QaQC _ QaQc^ 

the matrix elements c°''^d are in upper triangular form with zeroes in the diagonal when expressed in terms of the 
basis {K'^, Q", Q"^, L"^}- Therefore, the assumptions of Lemma [1] are satisfied and we obtain a well posed IBVP. 

C. Einstein's equations in harmonic coordinates 

As a second application of our theorem we consider Einstein's field equations in (generalized) harmonic coordinates. 
For this, we follow [l^, lla| and choose a fixed background metric g^^i, on M ~ [0,T] x S with the property that each 
time-slice St = {t} x S is space-like and the boundary T = [0,T] x 91] is time-like with respect to g^i,. We impose 
the following gauge condition on the dynamical metric gab, 

C' := g^' (r^t - r\b) -H' = 0. (10) 

Here, H'^ is a given vector field on AI and V^ab and Vab are the Christoffel symbols corresponding to the dynamical 
and background metrics, respectively. In the particular case where H'^ = and where the background metric is the 

Minkowski metric in standard Cartesian coordinates, V^ab vanishes, and the condition C^ — reduces to the usual 
condition for harmonic coordinates Dx'^ = for ^ = t,x, y, z. However, the advantage of the condition (jlOp is that it 
maintains the covariance of the theory since C^ is the difference between the two Christoffel symbols, 

C'ab = r^ab - T'ab = ^ff"" {"^ahd + ^bhad " ^dhab) , (11) 

where hab = gab — gab denotes the difference between the dynamical and the background metric. 
With the condition (fTU|) . Einstein's field equations are equivalent to the wave system 



g'^yAhab = 2g,fg'''C\cC^bd + ^ C^ dia9b)eC' cf g""^ - 2 g^"" R\d{a9b)e 

+ WttG (Xab ~ ^gabg'^'T.d] + 2V(aHb) , (12) 

°a , ^ 

where R bed denotes the curvature tensor with respect to (^^j^. Tab the stress-energy tensor and G denotes Newton's 
constant. Solutions of this equation which are smooth enough imply that the constraint variable Ca satisfies 

S^'^VcVdCa = -Ra^Cb - 167rGV^Ta6 . (13) 

Therefore, the imposition of the boundary condition Ca\q- = implies that any smooth enough solution of (|12p with 
initial data satisfying 

Ca|s„ = 0, n''VaCb\^, = 0, 

satisfies the constraint Ca = on M provided the stress-energy tensor is divergence free, V^Tab — 0. 

In order to formulate asymptotically nonreflecting boundary conditions we first construct an adapted local null 
tetrad {K°-,L°- ,Q°-,Q°-} and a radial function r as in the electromagnetic case. Notice that here these quantities 
are defined with respect to the dynamical metric gab and not the background metric g^j^. The boundary conditions 
derived in section HV CI are the following: 



LK'^K^'K^ir'^hbc) 

LK''K''LAir^hbc) + ^thbc 

^K-K'Q^ir^hbc) 
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= -QKK 


T 


- ^IQQ 


= 


= -qKQ 



(14) 
(15) 
(16) 



K'^Q^'Q^^ahbc - Q'^Q^'K^Vahb, 



T 



-QQQ > (17) 



(i^'^g^'O^ + L'^K^'K'^ - Q'^K^'Q' - Q'^K'^Q') Vahc = -2K''Ha\r , (18) 



r 



(if "L V + i"/^ V - Q"^'i" + 0"Q V) Va/^bc ^ -2Q''Ha\r , (19) 

(if"L''L^ + L'^Q''Q^-Q''Q''L"-Q"Q''L")4/j&c| = -2L''iJ,|^ , (20) 

where qxK and ggg are real-valued given functions on T and qkq and ggg are complex-valued given functions on T. 
The first three equations (jl4p .(fT5 |) . (fT6|) are related to the gauge freedom, the condition (fT7|) controls the gravitational 
radiation while the remaining conditions (fT51) . (|19p . (P0|) enforce the constraint Cq = on the boundary. The evolution 
equation (|12p has the form JT]) where E is the vector bundle of symmetric, covariant tensor fields on M and the 
boundary conditions (|14H20|) have the form Q where a = 1 and c°'^'^de is in upper triangular form when expressed in 
terms of the basis {K^K", K'^^ L''\ K'^^ Q''\ Q^ Q\ Q^^ Q''\ L^^ Q''\ L^ L"} . 

II. THE WAVE EQUATION ON A CURVED BACKGROUND 

In this section we prove Theorem [1] for the case of a single wave equation 

g'^^V.Vb^ = S (21) 

on M — [0,T] X S. For simplicity, we also assume that gab and 5 are independent of 0. In this case, it is convenient 
to choose Vq to be the Levi-Civita connection with respect to gab- The IBVP consists in finding solutions of ((2T|) 
subject to the initial conditions 

</.|s„=0o, n^V60|j.^=^o, (22) 

where 00 and ttq are given functions on Sq, and the boundary conditions 

[T''Vb<i) + aN''Vb<f\r^G, (23) 

where G is a given function on T. Here, n^ and N^ denote the future-directed unit vector field to the time-slices Et 
and the outward unit normal vector field to T, respectively, T^ is an arbitrary future-directed time-like vector field 
which is tangent to the boundary surface T and a is a strictly positive function on T . Without loss of generality, we 
assume that T" is normalized such that gabT°'T^ = — 1. Furthermore, by redefining cj) and S if necessary, we may also 
assume that the boundary data G vanishes identically. 

In order to show well posedness for this problem, we use a geometric reduction to a first order symmetric hyperbolic 
system with maximal dissipative boundary conditions [3,|j,[i3|- First, introducing the variables Va = Va(t>: the wave 
equation can be rewritten as the first order system 

Va0 = Va , (24) 

g'^'VaVb = S, (25) 

^aVb - VbVa - 0. (26) 

Next, we specify any future-directed time- like vector field m" and contract the first and the last equation with it. This 
yields the evolution system 

A0 = u^Va = n, (27) 

g'^'VaVb = S, (28) 

£uVb - VfcH, (29) 

where £u denotes the Lie derivative with respect to u°. This system is subject to the initial and boundary conditions 

4>\^o=(l^0, ^''^b\^^ = T^O , LoVb^ ioVb(pO , (30) 

[T''Vb + aN''Vb]^ = 0, (31) 

where to : Sq -^ M is the inclusion map, and subject to the constraint Ca = 0, where the constraint variable Ca is 
defined as Ca — Va — Va0. The evolution equations (|27p and (P5)l imply that Ca is Lie-dragged by the time evolution 
vector field u", 

£u^a = U. 



In the following, we assume that m° is pointing away from the domain at the boundary. This implies that a solution 
of (|27l28l29p with constraint-satisfying initial data automatically satisfies the constraints everywhere on M, and no 
extra boundary conditions are needed in order to ensure that the constraint Ca — propagates. 

Still, there is a huge freedom in choosing the evolution vector field u°; different choices lead to first order evolution 
systems (|27l28l29p which are incquivalcnt to each other if the solution is off the constraint surface Ca — 0- In this 
work we exploit this freedom in order to obtain energy estimates which allow for an appropriate control of the fields 
not only in the bulk but also on the boundary of the domain (see the estimate p6p below). In order to analyze this, 
following p/Tj we rewrite the evolution system (|28l29p in the form 

A^c^aV' = -u"(Val4 - VbK) + UbS/aV" = UbS, 

where the symbol is given by A'^bc = ~'^°5fcc + 2<5°(;,Mc). Since A'^bc is symmetric in be and since UaA°'bc = ~'^a'^°5fcc + 
2ubUc is positive definite, the evolution system is symmetric hyperbolic. In particular, the evolution equations imply 
that 

ya{A\,V'vn = {VaA\c)V'V^ + 2{ubV')S. 

Integrating both sides of this equation over the manifold M — [0, T] x E and using Gauss' theorem, one obtains^ 

naA\cV'V'^ = / naA\cV'V'^ + J NaA^cV'V^ ~ j [{V aA^cW^V^ + 2{ubV^)S] . (32) 

Et So T M 

The following two conditions (see 3) guarantee that the IBVP (|27l28l29l30l3ip is well posed: 

(i) naA°'bc is positive definite. 

(ii) For each p £ T, the subspace Af-{p) C TpM consisting of the vectors V''{p) satisfying the boundary condition 
(l?T|) at p is maximal non-positive. This means that NaA'^bc{p)V^ {p)V'^{p) < for all V^{p) € A/'-(p) and that 
A/L (p) does not posses a proper extension with this property. 

For the following, we choose the time evolution vector field u°- such that w" is everywhere future-directed and 
time-like on M and such that w° lies in the plane spanned by T" and iV° at each point of the boundary, more 
specifically, 

with < (5 < 1 a function on T. The following two lemmas imply the satisfaction of the conditions (i) and (ii) for an 
appropriate choice of S. 

Lemma 2 naA°'bc{p) is positive definite for all p £ AI . 

Proof. Let hab = 9ab + rLaUb be the induced metric on Sj and expand Ua = iJ-{na + Ua), where /i — ~n°'Ua. Since w" 
is future-directed and time- like, /i > and u°'Ua < 1. Therefore, 

naA^bc = M {hbc + ribUc + 2n(bUc)) 
is positive definite. r—i 

Lemma 3 Let Q < 5 < a(l -I- a^)^^ . Then, the boundary spaces M-(p) are maximal non-positive for all p (£ T. 
Proof, (cf. appendix B in Ref. 3) Fix a point p eT, and let V'' G TpM. We have 

NaA\cV''V^ = [6TbT, + 6NbNc + 2Tf^bNc)-SHbc]V''V^ 

= -5 [{T^'Vbf + (7VVb)2 + HbcV'V^] + 2 [S{T'Vbf + SiN'Vb)^ + (TV,)(7V-K)] , 



^ Notice that since n" is future directed, its flow increases t; hence in coordinates (t, x') where t parametrizes [0, T] and x' are local 
coordinates on S, we have n* > and rit < 0. 



where Hbc = Qbc + T^Tc — N^Nc is the induced metric on the orthogonal complement of the plane spanned by T^ and 
N^ . Eliminating the terms {T^Vb) in the second square bracket on the right-hand side using the boundary condition 
(ISTIl we obtain 



NaA^cV^V = ~5 [[T^Vbf + [N'^Vbf + i/fcc^V^] + 2 [5{a^ + l)-a\ [N'^Vbf 



(33) 



The last term on the right-hand side is non-positive by the assumption of the lemma. Therefore, NaA°'bc is negative- 
definite on the subspace of vectors V^ satisfying the boundary condition. Finally, we observe that M- (p) is maximal 
since its dimension is d = diuiTpM — 1 while the symmetric bilinear form NaA°'bc has signature (1, d). t—\ 



If we relax the assumption of homogeneous boundary data and replace the condition (|31|1 by the condition 

[T''Vb + aN''Vb]^^G, 



(34) 



we obtain, instead of ll33|) . 

NaA\cV^V'' = -5 [(TVb)^ + (iVVf,)2 + HbcV'V] + 2 [.^(a^ + 1) ~ a] (N^Vb)'^ + 2(1 - 25a)iN^Vb)G + 25G^ 
Let < p < 1 and set 6 — {1 — p)a(l + a^)^^. Then, we have (cf. appendix B in Rcf. [9]) 

(1 - 2(5a)2 



NaA\cV''V' < -6 [(TVb)2 + {N''Vb)^ + HbcV'V'] + 



25 



2ap 



G' 



(35) 



Using this and the positivity of naA°'bc in the identity 



/ UaA^cV^V < I UaA^cV^V - Ci f naA\cV''V' + C2 I G 



we obtain the estimate 
t 
C3 



St 



is 



naA\aV''V' + / S' 



ds 



for all < t < r, where Ci , C2 and C3 are strictly positive constants which are independent of V^, and Tt := [0, t] x 3E. 

t 
Applying Gronwall's lemma'^ to the function y{t) '■= J J UaA'^bcV^V'^ds we obtain from this 

s.. 

Lemma 4 Let T > 0. There is a constant C = G{T) > 1 such that all smooth enough solutions to the IBVP 
\28\29\3(]y34^ satisfy the inequality 



naA^cV'V' + / naA^cV'V' < C 

St Tt 

for allO<t<T, where % := [0, t] x ^E. 



naA\cV^V' + G^ 

ISo Tt 




(36) 



Since any solution of this problem also satisfies u°'Ca = u°-Va — £u<l> = 0, £uGa = and IqCq = Lo{Va — ^ a4>) = Oi 
and since w" points outward from the domain at T, the constraint Ca = is satisfied everywhere on M . From this 
and the previous lemma, we have established: 



Theorem 2 The second order problem Ii21l22l23\) is strongly well posed: given smooth initial and boundary data <j)o, 
TTo and G satisfying the usual compatibility conditions at 9So, there exists a unique smooth solution satisfying the 
estimate \3b]) with V"" replaced by V°(/). 

Remark 1 The important feature of the estimate iSb]) is the second term on the left-hand side which yields a 1? 
boundary estimate for the gradient of </>. This estimate is obtained by choosing the time evolution vector field w° in 
such a way that the boundary matrix NaA°'bc is negative definite on the subspace of vectors satisfying the boundary 
conditions. As we will see (Lemma\^ in the next section), this property is important for systems of wave equations 
since it allows the coupling of the boundary conditions through small enough terms involving first derivatives of the 
fields. If, on the other hand, u"" is chosen to be tangent to the boundary, the boundary matrix has a nontrivial kernel 
and one does not obtain an estimate for the gradient of 4> on the boundary from the first order system. However, this 
does not affect the strong well posedness of the second order system which is independent of u°' . 



See, for instance, Lemma 3.1.1 in Ref. [IS 



As an example, consider the wave equation on the half-plane S = IR+ x R^ with the flat metric g = —dt^ + dx^ 
dip + dz^. In this case, we have 

p 



with {pyf + (/3^)2 < 1 and p := ^1 - (/3!^)2 - {I3'Y , and the boundary condition ([23]) reduces to 

[0t + pac^^ - py^y - /3^0.],^o = pG, (37) 

where (^^ := 9f0 etc. Choosing li" = p{T'^ + ^A^"^) with < ^ < 1 the energy norm for this problem reads 



CX3 OO OO 



naA\cV''V^ = / / / [^^t+ ^1 + ^1 + ^1 + '^<f>t {5pc^. + py(j>y + /3"0.)] dy dz dx. 

St — oo — oo 

This is similar to the norm we used in Ref. jB] for obtaining an a priori energy estimates for the second order wave 
equation with boundary condition (jST]) . 



III. SYSTEMS OF WAVE EQUATIONS AND PROOF OF MAIN THEOREM 

In order to show that the system (|1I2I3P yields a well posed IBVP, we follow the arguments given in Sect. |TT] 
and reduce it to a first order symmetric hyperbolic system with maximal dissipative boundary conditions. Let 
Va'^ := Va<&'^, and let u°(p, $) denote a future-directed time-like vector field on M such that 

with < (5 < 1 a function on T to be determined. Then ([T]) can be rewritten as the first order evolution system 

u"'Va<^^ = W'Va^, (38) 

g''\<^)^aVb^ = S^{<^,V), (39) 

u'^i^aVb^-V.Va^) = U-R^Bab'^''. (40) 

where R Bab denotes the curvature belonging to the connection Vq. At this point, we stress that the connection Vq 
is a fixed background connection on the vector bundle E, and not the Levi-Civita connection belonging to the metric 
gab{^)^ so that R^ Bab does not depend on <I> nor its derivatives. The system (|38l39l40p is subject to the constraint 
Cb"^ = 0, where Cb^ := Vh$"* — Vb^ ■ Eqs. (|38I40|) imply that the constraint variable Cb^ is Lie-dragged by u": 

£„C/ = u-^VaC/ + (Vfcu'^)a^ = 0. 

Therefore, any smooth enough solution of the first order problem p8l39|40p belonging to initial data with C\/^ — 
satisfies the constraint Cf, = everywhere it is defined. The initial condition is 

^^\^^^^^, nV,^|^^=^^ tlV,^ = Ll^^^t. (41) 

and the boundary condition ([3]) reads 

[rVb + aTVVb]^ = e^B Va''\^ + d^B '^''\r + G^- (42) 

In order to analyze the well posedness of the first order IBVP p8|391 I40|41|42]l we first linearize the system by 
replacing the coefficients 5a6($), S''^($, V$), T''($), Ar''($), a($), c°-^b{'^): d^si^) by smooth functions 5a6, S^, T^ 
7V^, a, d^^Bi d"^B, respectively. Local in time well posedness for the original quasilinear system follows by iteration 
from the well posedness result for the linear system with enough differentiability^. Next, we use a partition of unity in 
order to localize the problem. With this, it is sufficient to consider a local trivialization ip : [/xR^ i— s- 7r~^([/) of _E such 
that U C M is compact and contains a portion U of the boundary T. Let e > 0. According to the assumption there 



** See, for instance, [lOl. llSH . 



exists a smooth map J^ : U ^ GL{N,M.),p i^ {J^{p)) such that the transformed matrix coefficients c° := J^c°'J~^ 
satisfy the condition @ for aU vector-valued one-forms Va onU. Setting /lyis(e) :— {j]'hJ^)AB — hcoiJe)'" a{Je)^ b, 
we can reformulate this condition by stating that 

hAB{eV^ci^)c''^ni'^)Va''V,^ < ehMe)e''\<^)Va^V,'' , (43) 

for all vector- valued one- forms V^ on U. The system (I38I39I40I) can be written in the form 

-AhAB{e)u-Va \ f '^'^ \ - s(<P V) f44^ 

where A > is to be determined, A"'bc — —u°-gbc + 25"(6Wc) and 



-KhAB{e)u''Va'' 
hAB{e)R^Cab^'^'u- + hABie)ubS^{^, V) 






Let B{na\ (^, W), ($, V)) denote the bilinear form belonging to the principal symbol of (f44|) . that is, for an arbitrary 
one-form Wa on M define 

B{Wa; (*, W), ($, V)) := -Ku'^WahABie)-^^^'' + hAB{e)WaA\cW''^V^''. 

We have 

Lemma 5 Let A > 0. Then, B(na', (^, W), ($, V)) is symmetric in (^, W), ($, V) and positive definite for Wa — Ua 
and Wa — ria- Therefore, the system |^^[) is symmetric hyperbolic. 



Proof. The symmetry property follows immediately from the symmetry of /i^s(£) and the symmetry of y^"t,c in be. In 
order to check the positivity statements, let Wa — Ua, 7 := \J~-vFu^ and Ua '■= l^^Ua. Since A°'bcUa — 7^ [dbc + Z-UfcUc], 
we find 



BK; ($, V), ($, V)) = 72 [A/i^s(£)$^$^ + {gab + 2uaUb)hAB{e)VaAV''''] 

which is manifestly positive definite. The proof that B{na] ($, V), ($, V)) is positive definite is similar to the proof 
of Lemma [51 r—i 

As in the previous section we obtain well posedness of the linearized system provided we can show that each 
boundary space 

AL(p) := {($,y) e R^ X m('^+i)^ : [t\p) + a(p)N\p)\ 1// = c'^^ BipWa"" + d^ b{p)<^''}, peU, 

is maximal non-positive with respect to B{Na\ ($, V), ($, V)). This is the statement of the next lemma. 

Lemma 6 Set 5 :— a{\ + a^)^^ /2 and n :— 2[26 + (1 — 25aY / a^ . Choose e > small enough such that ke < S and 
A > large enough such that 2KhAB{£)d^cd^ D^'"'^^ < (5A/iAB(e)$'^*^ for all $ e R^. Then, the boundary space 
J\f-{p) is maximal non-positive for all p ^lA. 

Proof. Let p ^U. We have, as in the proof of Lemma [3l 

B(7Va;($,1/), ($,!/)) = -Au''NahAB{s)<^^<^''+hAB{£)NaA\eV''^V' 

= -ShABie) [(T°T'' + N'^N'' + H''^)Va^Vb^ -f A$^$^] 

-f 2[5T''T^ + 5N''N^ + T''N^]hAB{£)Va'^Vb^. (45) 

Let ($^,Va^) € AA(p). Then, T'^K^ = -aN'^Va'^ + G^ with G^ := c'^'^BVa^ + d^B^^, and we may use this 
equation in order to eliminate the terms [T°-Va^) in the second bracket on the right-hand side of (|45)) . This yields 



B{Na; ($, V), ($, V)) < -5hAB{e) [(T^T'' + N'^N^ + H''^)Va^Vb^ + A*-^*^] 



26 



(1 - 2Sa) 



21 



/iAB(£)G^G^, 



where we have set 6 :~ a{l + a^) ^ /2 and used the boundary estimate (fSS]) with p = 1/2. Now, 



/iAB(£)G^G^ < 2hAB{£V^cVa''c''''DVb'' + 2hABi£)d^C^''d''D^'' 

< 2ehAB{£)e'"'Va^Vb'' + 2hAB{£)d^cd''D<^''<^'', (46) 
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where we have used the estimate (gS]) in the last step. Recahing that e"'' = g"'' + 2T''T'' = T'^T'' + N'^N'' + H"'' and 
the definition of k in the assumption of the lemma we find 

B{Na;{^,V),{<P,V)) < -ShABie) [e^'Va^n'' + Al'"'^''] + >i[ehABie)e'''Va^Vb'' + hAB{e)d''cd'' d^''^''] ■ 

The non-positivity oi Af-{p) now follows from the assumptions on e and A. Finally, we observe that an element in 
M-(p) is characterized by N conditions in a (d + 2)A^-dimensional space which implies that dm\M-{p) > {d-\- 1)N. 
On the other hand, from Eq. (|45| we see that the signature of B{Na\ ., .) is given by [N^ {d + l)N). Therefore, 
dim J\f_{p) = (d + l)N and the maximality oi J\f_{p) follows. r— i 



IV. BOUNDARY CONDITIONS FOR ISOLATED SYSTEMS 

We consider here boundary conditions for an isolated system emitting radiation. If, for computational purposes, 
the evolution domain of such a system has a finite (artificial) boundary, some artificial boundary condition must 
be imposed. If one knew the correct boundary data for the analytic problem, then in principle one could use any 
boundary condition corresponding to a well posed IBVP. However, the determination of the correct boundary data 
is in general a global problem, in which the boundary data must be determined by extending the solution to infinity 
either by matching to an exterior (linearized or nonlinear) solution obtained by some other means. The matching 
approach has been reviewed elsewhere [1!||. Here we consider an alternative approach in which homogeneous boundary 
data can be assigned in such a way that the accuracy of the boundary condition becomes exact in the limit that the 
boundary is extended to infinity. (Such boundary conditions would also be beneficial to the matching approach 
because the corresponding boundary data would be small so that numerical or other error would also have a small 
effect.) Artificial boundary conditions for an isolated radiating system for which homogeneous data is approximately 
valid are commonly called absorbing boundary conditions (see e.g. |20l. l2ll. |22| . l23l. l24l 1 25J ) . or nonreflecting boundary 
conditions (see e.g. [2^, [23, ^^) or radiation boundary conditions (see e.g. Q). Such boundary conditions are 
advantageous for computational use. However, local artificial boundary condition are not perfectly nonrefiecting in 
general. Here, to be more precise, we consider nonreflecting boundary conditions in the sense of boundary conditions 
for a well posed problem for which homogeneous data produces no spurious reflection in the limit that the boundary 
approaches an infinite sphere. The extensive literature on improved versions of nonreflecting boundary conditions 
involves higher order and nonlocal methods. Our interest here is to investigate the optimal choice of local first 
order homogeneous boundary conditions on a spherical boundary for the constrained Maxwell and linearized Einstein 
problems expressed in terms of the gauge dependent variables A'^ and j^'^ . See [30l l3ll . |32 | for the construction of 
higher-order and higher-accurate boundary conditions for Einstein's equations. 

We base our discussion on waves from an isolated system satisfying a system of fiat space wave equations. We 
use Greek indices to denote standard inertial coordinates x'^ = (t, x, y, z) in which the components of the Minkowski 
metric rj^''' are diag(— 1, 1, 1, 1). In the case of a scalar field $, we thus consider the wave equation 

v'"^do.dp<^ = {-df + dl + d^ + 92) $ = 5, 

where the source S has compact support. Outside the source, we assume that the solution has the form 

lit rv*^ jytfj 

where (r, 6, (f) are standard spherical coordinates and /, g and h and their derivatives are smooth bounded functions. 
These assumptions determine the exterior retarded field of a system emitting outgoing radiation. The simplest case 
is the monopole radiation 

r 
which satisfies [dt -t- dr){r^) — 0. This motivates the use of a Sommerfeld condition 

-{dt+dr){r^)\R=q{t,R,e,cj,) 
r 

on a finite boundary r — R. 
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The resulting Sommerfeld boundary data q in the general case (|47)) falls off as l/i?"^, so that a homogeneous 
Sommerfeld condition introduces an error which is vanishingly small for increasing R. As an example, for the dipole 
solution 

'^Dipole — Oz = ' 

r 
we have 

f{t — r) cos 9 



q 



i?3 



A homogeneous Sommerfeld condition at r = _R would lead to a solution ^oipoie containing a reflected ingoing wave. 
For large R, 

~ F{t + r-2R)cos9 

^Dipole ^ ^Dipole + K , 

r 

where dtf{t) ~ F{t) and the reflection coefficient has asymptotic behavior k = 0{1/R^). More precisely, the Fourier 
mode 

/ piuj{t-r) iu](t+r-2R) 
^Dipolei^) = d^ I h Kuj 

\ r r 

satisfies the homogeneous boundary condition {dt + dr){r^Dipoie{'^)\R — with reflection coefficient 

1 1 



(48) 



" 2cj^R^+2iu;R-l 2lj^R^' 

Note that (59) and (62) satisfy 

K ~ qR. (49) 

In the case of a system of equations k will have N components corresponding to the number of modes generated in 
the reflected wave. The boundary conditions lead to a system of simultaneous equations relating k to the components 
of the Sommerfeld data q. If these equations are nondegenerate then (|49| continues to hold. However, degeneracies 
could conceivably lead to weaker asymptotic falloff of k. (It would be interesting to determine whether such cases 
exist.) In any case, (|49|l gives the optimum allowable behavior of the reflection coefficients so that the asymptotic 
behavior of the Sommerfeld data q is a good indicator of the quality of the boundary condition. This forms the basis 
of our investigation of the Maxwell and linearized Einstein equations with a spherical boundary in sections IIVBI and 
HVCl 

A. A plane boundary 

The key ideas in the above example are that (i) the Sommerfeld condition is only satisfied exactly by waves traveling 
in the radial direction and (ii) in the asymptotic limit r -^ oo all waves from an isolated system propagate in the 
radial direction. This allows us to reformulate our discussion of the Sommerfeld condition by considering a wave $ 
propagating in the domain a; < 0, which is incident on a plane boundary at a; = with the boundary condition 

where K^da — dt + dx is the characteristic direction determined by the outward normal to the boundary dx and the 
time direction dt- This homogeneous condition is satisfied for plane waves $ = G{t + k^x + kyy + fc^z) incident on 
the boundary only for the single case {kx,ky,kz) = (1,0,0), i.e. a plane wave propagating in the outgoing normal 
direction. Plane waves in the normal direction pass through the boundary whereas plane waves incident in other 
directions on the boundary give rise to a reflected wave. We will take advantage of this simplification of the plane 
wave case in discussing boundary conditions for electromagnetic and gravitational waves. The results then suggest how 
to formulate boundary conditions for an isolated electromagnetic or gravitational system with a spherical boundary 
of radius R, where in the limit R —f oo all radiation is incident normally. 
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For the electromagnetic case, we describe the field by means of a vector potential A^^ satisfying the Lorentz gauge 
condition. Maxwell's equations in a flat spacetime with Minkowski metric r]^'^ then reduce to the wave equations 

?7"^9„5/3A^ ^ 

subject to the constraint 

C := df.A" = 

introduced by the Lorentz gauge condition. This constraint keeps us from requiring that each component of A^^ satisfy 
a homogeneous Sommerfeld condition, in contrast to the scalar example. The electromagnetic case also differs from 
the scalar case because of the remaining gauge freedom allowed by the Lorentz condition. 

An electromagnetic plane wave incident in the outgoing normal direction can be described by the real part of the 
vector potential 

A^ = F{t - x)Q^ + G{t ~ x)K^ , 

where F{t — x) is complex, Q^ = Y^ + iZ^^ is a complex null polarization vector, G{t — x) represents gauge freedom 
and K'^ = T^ + X'^, in terms of the orthonormal tetrad (T^, X'^, V^, Z^) aligned with the coordinate axes satisfying 

Vfj.i' = ^T^Ty + X^Xi, + Y^Yy + Z^Zy . 

In order to formulate a gauge invariant boundary condition we consider the corresponding electromagnetic field 
tensor 

i^p. = ^^,A, ~ dyA,, = ~F'{t - x){K^Qy - Q^K,). 

Here we adopt the notation duF{u) = F'{u). For this plane wave, all components of -F^^ satisfy 

K^'F^, = 0. 

However, this condition rules out the possibility of a static electric field oriented normal to the boundary. For the 
purpose of formulating a boundary condition which only restricts propagating waves it suffices to consider the weaker 
condition 

K^Q'F^y = 0. (50) 

In terms of the electric and magnetic field components tangential to the boundary, (j50p corresponds to the plane 
wave relations Eta„ • 'Btan — and |E4a„| = \'Btan\, with the corresponding Poynting vector in the outward normal 
direction. 

We can incorporate (|50|) into the following homogeneous Sommerfeld boundary conditions for the vector potential: 

K'^Kt'd^A, = 0, (51) 

Q'^Kf'd^A, = K'^Q^'d^A,. (52) 

The remaining boundary condition can be expressed in Sommerfeld form by rewriting the constraint as 



C=- {-L^K^ - K^L^ + Q^'Q^ + Q'^Q^) d^^A, = 0, (53) 



where L^ = T^ — X^ . Here [K^ ^ L^, Q^) form a null tetrad according to the conventions 

nt,v = -K(^L^) + (3(^(3^). (54) 

We assume throughout the following that the spin transformation freedom Q'' — > e*"^'' has been restricted according 
to K'^d^a — 0. The Sommerfeld boundary conditions (fSTjl . ([52|l and ([53|l have the required hierarchical, upper 
triangular form for a well posed IBVP, see Lemma [TJ 

For the purpose of extending this approach to the gravitational case, we write the linearized Einstein vacuum 
equations in the form 

V'^^dadp^'' = (55) 
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subject to the harmonic constraints 

C^ ■■= -d^^l^" = 0. (56) 

Here, to hnearized accuracy, we set y/—gg^^ = rj^'^ + 7'"' so that ^^v represents the perturbation of the densitized 
metric. (Indices of hnearized objects are raised and lowered with the Minkowski metric.) 

A plane gravitational wave incident on the boundary in the outgoing normal direction is given by 

7^^^ = F{t - x)Q^'Q'' + K^^^C\t -x)- ^n^^'^K'^^^it - x), 

where ^'^ {t — x) describes the gauge freedom. 

In order to formulate a boundary condition with gauge invariant meaning analogous to (j50p in the Maxwell case, 
we consider the linearized curvature tensor 

- 2i?p^pCT = dpd^jf^a - dad^jf^p - dpdf^-fycr + d^Opj^p - -{ri^iaOpdiyj - rippdad^-f - ij^adpdp'y + rj^pdadpj) 

= 4F"{t-x)K[pQ,]Q[pK^^ . (57) 

Plane wave boundary conditions on the curvature tensor could be imposed by requiring that the Newman-Penrose 
component \E'o = K^Q'^Q^K'' Rpypcr vanish on the boundary. (See [2] for a discussion of the appropriateness of this 
boundary condition.) However, this requirement involves second derivatives in the normal direction when expressed 
in terms of ^p^. Instead, we require \1/ := K^Q'^Q^T"' Rp^pa^ = on the boundary. The condition ^ = is equivalent 
to ^0 = if the Ricci component RpvQ^Q'^ = 0, e.g. if the vacuum Einstein equations are satisfied. 
A straightforward calculation leads to 

-2* = K^Q''QPT%dpd,^p,-d^d,^pp-dpdp^,^ + d„dp^,p) + ]^Q''QPd,dp^ 

- K^Q'Q'T^i-d^d^^pp - dpdp^,, + d^dp^^p) + \^(K^^K- + Q^Q'')Q'QPd,dp-ip„ 

= Q^d, (\{K^K^ + Q''Q")Q''dp-ip, - K^QPT'^dp^^p 



+ T^d, ( - K^^Q^Q^duipp + K'^Q^'Q^dp^^p] . (58) 

Thus, besides containing no second derivatives normal to the boundary, the condition ^ = can be reduced to two 
first order conditions by factoring out the Q'^di, and T'^d^r derivatives in ([55)1 which are tangential to the boundary. 
There are many ways this can be done. In order to obtain first order conditions which fit into a hierarchy of Sommerfeld 
conditions, we modify (|58)) according to the steps 

i(if^i^- + Q^QnQ'dpip, - ^RpQ^L-dpjp, ~ i. 



2* = Q'd, -{K^^K- + Q^Q^Q^dp^p, - -K'Q^^L'^dpjp^ ~ T:K^Q''K''dpj,p 



T^d, - K>'Q''QPd,-ipp + K^'Q^'Q^dp^^p (59) 



= T^Q'du {Kf^K^QP + Qf'K^LP - Q^Q''QP)dp^p^ - IQ^Cp - K^QPR^Opj^p 



1 

-( 

2 

+ T^d, I - K^Q-'QPd^jpp + KPQ-'QPdp^^p ) (60) 



= -Q'^d, ( (KPR^QP - Q^Q^QP)dp'ip, - 2KPQPK^dp^^p - 2Q^Cp 



1 

-( 

2 

+ T^dAKPQ''QPdp^,p\. (61) 

Thus since the derivatives Q'^di, and T'^di, are tangential to the boundary, we can enforce ^ = on the boundary 
through the first order boundary conditions 

Q^QPRPdp-f^p = 0, (62) 

K"Q^KPdp-/a,p - ^K^K'^Q'^dpjo.p + iQ'^Q^QPdp-io.p = 0. (63) 
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These two boundary conditions can then be included in a hierarchical set of Sommerfeld boundary conditions, ac- 
cording to the example 

K^K^K^'d^-i^p = 0, (64) 

Q'^Q^K^d^lcp = 0, (65) 

Q'^Q^K^d^^^p = 0, (66) 

K''QPK^^^^^^p - ^K''K''Q^d,,-f^0 + ^Q''Q>'Q>'d,nc.t3 = 0. (67) 

The constraints Cp — 0, which determine the remaining boundary conditions, can be cast in the Sommerfeld form 

which can also be incorporated into the hierarchy. 

However, there are many alternative possibilities to ([M)l - (|S7|) which preserve the hierarchical Sommerfeld structure 
and lead to a well posed IBVP. In the absence of a clear geometric approach, we next examine the boundary conditions 
appropriate to an isolated system by considering the resulting reflection off a spherical boundary. 

B. Application to Max'well fields ■with a spherical boundary 

In the case of a general retarded solution for a massless scalar wave equation, we found that a Sommerfeld boundary 
condition on a spherical boundary of radius R required data q = 0(1/ R^). Homogeneous Sommerfeld data gave rise 
to an ingoing wave with reflection coefhcient k = 0(l/i?^), as in (P5|) . This is the best that can be achieved with 
a local first order homogeneous boundary condition on a spherical boundary. We now investigate the corresponding 
result for the constrained Maxwell equations expressed in terms of a vector potential A^. 

In doing so, we associate spherical coordinates (r, x ), x = {9, 0), in a standard way with the Cartesian coordinates 
x* = (x, y, z), e.g. z = r cos 6. As in (|54p we introduce a null tetrad {K^, L^, Q^) adapted to the boundary, where now 
K^dp — dt + dr, Lf^dp — dt — dr and we fix the spin-rotation freedom in the complex null vector vector Q^ — (0, Q^) 
by setting 

^' ^''\q^. (68) 



where 



dx^ 



Q^=(Q«,Q^) ^1(1, -J 



r \ smt 



We describe outgoing waves in terms of the retarded time u = t — r. 

In order to investigate the vector potential describing the exterior radiation field emitted by an isolated system we 
introduce a Hertz potential with the symmetry 



Then the vector potential 



4 



A^' = d^H'"' 



satisfies the Lorentz gauge condition and generates a solution of Maxwell's equations provided the Hertz potential 
satisfies the wave equation. The trace H represents pure gauge freedom. 

We consider outgoing dipole waves oriented with the z-axis. Other dipole waves can be generated by a rotation. 
Higher multipole waves can be generated by taking spatial derivatives. 

The choice H = Z"da—^, Tjl'^''] = Q gives rise to the dipole gauge wave 

/F"{u) F'{u)\ ^^^ f2F'{u) 3F{u)\ ^^ f F'{u) F{u)\ ^ 

Ap - — ^ + -\^ cosOKp + ^ + — ^ cos^S^r - -\^ + —^ Zp 
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with components 



'F'{u) 2F{uy 
K^'A^ = ( — V^ H ^ I cost 



-A„.|£;w+£M|„„,. ,69, 



In appendix |^ we give some useful formulae underlying the calculation leading to (j69p and the following results 
The choice H^^^ = {T^^Z" — ^My^iW. gjygg pjgg ^q g^ dipole electromagnetic wave 



with components 



A.-['^.mT,.^^-£^z, 



At'K^ = ^cos6 



A^'Q^=^-^sine. (70) 



,/(") 



The choice W^^ — [X^^Y" ~ y^^i^jLLJ. gives rise to a dipole electromagnetic wave with the dual polarization 

A.,,= 



f'{u) , f{u)\fyX^ xY^ 



with components 



\ ^ 



A^K^ = 0, 



Ai^r, ■('/'(«) f{u)\ . , , 



We wish to formulate boundary conditions which generalize the Sommerfeld hierarchy (|5ip and (|52p to a spherical 
boundary of radius _R in a way which minimizes reflection. By inspection of (|69p . (|70p and (|7ip . we consider the choice 

^K'^d^ir^K'^A,) = qK , (72) 

^K^d^irQ-'A,) - Q'^d^iK-'A,) = qq, (73) 

chosen to minimize the asymptotic behavior of the Sommerfeld data. As before, the constraint determines the 
remaining boundary condition as part of the Sommerfeld hierarchy. 
For the dipole gauge wave (|69l) . 

2F(m)cos6' 
IK = ^4 , 9Q = 0; 



for the dipole electromagnetic wave ((701) . 






QK^O, qg = 3 sm6^. 



and for the dual dipole electromagnetic wave ((7T|) 

Overall this implies qk — 0{1/R'^) and qq — 0{1/R^). We have checked that homogeneous Sommerfeld data leads 
to reflection coefficients with overall behavior n — 0{1/ E?) in accordance with (|1^ . 
Note that the relations (|A1[) and (jA8[) allow us to express (|72p and ((73)) in the form 

^K'^K^^df^ir^A,) = qK , (74) 

Q''K^^^^A, - K-Q'^d.A. = Qq, (75) 

which correspond to ([7]) and ((S]) when d^ is generalized to the connection Va in a curved space background. Here 
ni is equivalent to the gauge invariant condition 

Q^K^F^, = qg. (76) 
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C. Application to linearized gravitational fields ■with a spherical boundary 

The gravitational case is more complicated than the electromagnetic case because the geometry of the boundary is 
coupled with the boundary condition. Additionally, there are no gauge invariant quantities, analogous to (j76p in the 
electromagnetic case, on which to base first order boundary conditions. We begin with a discussion of how to adapt 
to a curved boundary the first order version of the ^ boundary condition given in Sect. IIV Al for a plane boundary. 

In the nonlinear treatment of a curved boundary with unit outer normal TV" we can decompose the metric according 
to 

gab = Tab + NaNb , 

where Tab is the metric intrinsic to the time-like boundary. Let Da denote the covariant derivative associated with 
Tab- The extrinsic curvature of the boundary is 

Nab = Ta'^VcNb . 

We complete an orthonormal basis by setting 

Tab = -TaTb + Q(aQb) 

in terms of a time-like vector T" and complex null vector Q" tangent to the boundary. 

We decompose * := K"-Q^Q''T'^Rabcd = *t + *w and the Weyl component *o = K'' Q'' Q" K"^ Rabcd = *t + *7V + 
2^TN, where JsT" = T" + iV and 

*^ = T'^Q^'Q^T'^Rabcd , (77) 

^^ = N'^Q'^Q^T'^Rabcd , (78) 

*TJV = T^-Q'^Q^N'-Rabcd ■ (79) 

When the vacuum Einstein equations are satisfied the Ricmann curvature tensor may be replaced by the Weyl tensor 
whose symmetry implies ^tn = 0. Therefore, in this case, ^ = implies the vanishing of the Newman-Penrose Weyl 
component ^o = 0. 

A short calculation gives the embedding formulae 

^N = Q^Q'T\DdNbc - DbN,d) 
and 

*j, = T'^Q^'Q^T'' {^''^Rabcd - NacNbd + NbcNad) , 

where ^^'Rabcd is the intrinsic curvature to the boundary, i.e. 

T'^Q'^Q'T'^'^^^Rabcd = Q'^Q^T^DdD, - D,Dd)Tb . 

(These are the embedding equations for the Cauchy problem corrected for the space-like character of the normal to 
the boundary.) 

We now apply these results to a spherical boundary r = i? in linearized theory off a Minkowski background, i.e 
dtiii — Tj^ii-'rehp^i, in standard inertial coordinates a;^, where e is the linearization parameter. We choose T^ — dfj,t+0(e) 
and Nfj_ = d^r + 0{e). Then D^jT^ = 0(e) and iV^^ = R~^Q^^ + 0{e), where Q^^ = Q(^lQu) is the metric of a 2-sphere 
of radius R. We choose the basis to satisfy T^^D^T^ = and T^^D^Q^ = 0, so that 

*T = T^'Q''Q''T^<-''^R^,p, + 0(e2) = T" D^iQ^QPOpT,) + ©(e^) 

and 

*Ar = T^D^iQ'^QPNp,) - QPOpiQ^T^N,,) + W{DpQ^)QPQ''T''N„, + ^Q^Q^DpT, + ©(e^). 

Thus the boundary conditions 

Q'Q'iNp. + DpT,) = 0, 

Q'TPNp, = 0, (80) 
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imply to linearized accuracy that 



1 



* = -BQ'Q'DpT, ■ (81) 



This gives a geometric formulation of the first differential order version of the requirement that ^ — > in the asymptotic 
limit R -^ oo. However, vI'q — 0{1/R^) in an asymptotically flat space-time, whereas ([5T|) leads to ^ = 0{l/R^). 
This is an indication that the boundary conditions ((80|) might lead to more reflection than desirable. Can this be 
remedied by the introduction of, say, lower order terms in the boundary conditions? We investigate this question in 
the context of a well posed IBVP based upon the harmonic version of the linearized Einstein equations ([55)) and (|56p , 
where 7^" = -ht"" + ^•q'^'h. 

For this purpose, we now consider linearized outgoing waves in the harmonic gauge which arc incident on a spherical 
boundary. We model our discussion on the Maxwell case by using the gravitational analogue of a Hertz potential 
jjtiauf3 J3^^ [gj ^ which has the symmetries 

and satisfies the flat space wave equation 



d^d^H'"^"^ = 0. 



Then the densitized metric perturbation 



satisfies the linearized Einstein equations in the harmonic gauge. Outgoing waves can be generated from the potential 

r 

and its spatial derivatives. 

The incidence of such an outgoing wave on a boundary r ^ R leads to refiection, with the asymptotic falloff of 
the reflection coefficients depending upon the choice of boundary conditions. We limit our calculation of reflection 
coefficients to the case of outgoing quadrupole waves, which can be obtained from the Hertz potential 

j^tiavfi _ j^^auP JjW_ (82) 

r 

where K^°"^^ is a constant tensor. (All higher multipoles can be constructed by taking spatial derivatives.) K^^°"^^ has 
21 independent components. However, the choice K'^"'^^ — e^°"^l^ leads to 7'*" = so there are only 20 independent 
waves. These can be further reduced to pure gauge waves, corresponding to the trace terms in K^^°"^^ , e.g. A'^""'^ = 
j^ai/^/3M _ jyM'^^a/s leads to a monopole gauge wave. Linearized gravitational waves arise from the trace-free part 
of K^"-'^^ . There are ten independent quadrupole gravitational waves, corresponding to spherical harmonics with 
(£ = 2, —2 < TO < 2) in the two independent polarization states. The other ten independent potentials comprise two 
monopole gauge waves, three dipole gauge waves and five quadrupole gauge waves, for which the linearized Riemann 
tensor vanishes. It suffices to consider the following examples of waves with quadrupole dependence aligned with 
the z-axis. Other quadrupole waves can be obtained by rotation and have similar asymptotic behavior. Reflection 
coefficients from the other monopole and dipole gauge waves are smaller and provide no further useful information. 
The Hertz potential ([5^ gives rise to the perturbation 

r 
Appendix VK\ lists useful formula for the calculations underlying the following results. 

1. Quadrupole-monopole gauge wave 



The Hertz potential 



jjt^aup = fzf'ri°"'Z'^ + Z^rff'Z" - Z^'tj^'^Z" - Z'^ri'^'^Z' 



f{u) 
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gives rise to a combination monopole-quadrupole gauge wave with components 



1-fi = -2(^^ + ^jsin^0, 
K"Q^jc,i3 = -^^ sin 6* cos 61, (83) 



Here the sin 9 dependence of the spin- weight 2 component Q'^Q^^ap is a pure 2^20 spin-weighted spherical harmonic; 
the sin cos dependence of the spin-weight 1 component K^Q^jap is a pure 11^20 harmonic; and the remaining spin- 
weight components are mixtures of Yqq and 1^20 ■ 

2. Quadrupole gravitational wave 
The trace-free Hertz potential 



jjtiavp _ I cyM^a _ Z^T^^iX^Y^ — y^ X^) + {X^Y°' ~ Y^X°'){T'^ Z^ — Z'^T^) | 



(84) 



gives rise to a perturbation with 7 = and components 

Q"Q''7a/3 = 2zsin 



«n.. . _o...2.^/"H , /'(-) 



Q^Q^lcff = 0, (85) 

i^"Q^„, = ^cos0sin-^'^'(") ■ ^^^"^ 



V *" 

which have spin-weighted ([ — 2, m — Q dependence. 

3. Dual quadrupole gravitational wave 
The trace-free Hertz potential 

jjt^c^^p = ({T^'Z°' - Z^'T°'){T'' Z^ - Z^T^) - [X^'Y" - Y^'X°'){X''Y^ - Y" X^) + -{ti'^'tj"^ - t?^'^??''") j 
obtained from the dual of (|84p . gives gives rise to a perturbation with 7 = and components 

Q-g^., = 2sin^«^/>)^/>)-.M 



/H 



rg^., ^4(cos^.4)(/M + M), (86) 



if"if'^7„, = 2(cos2 0-i)M 

3 r'^ 



which have spin-weighted £ = 2, m = dependence. 
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4- Sommerf eld-type boundary conditions 

Sommerfeld boundary conditions consistent with a well posed harmonic IBVP have wide freedom regarding (i) 
partial derivative terms consistent with the hierarchical upper triangular structure of the boundary condition and (ii) 
lower differential order terms. Here we consider three choices of of boundary conditions and compare their reflection 
coefficients. One basic idea common to these choices has already be used in the scalar and Maxwell cases, i.e by 
inspecting the asymptotic behavior of the waves ([55]) . ([55]) and ([55)1 we use the property K°'daf{u) = to introduce 
the appropriate powers of r that lead to the smallest asymptotic behavior in the resulting Sommerfeld data. 

Our first choice of boundary conditions is the mathematically simplest choice 

^K'^KPK^'^^,{r''^^|i) = qKK , (87) 

WQ^K^^d^in^p) = qQQ , (88) 

^Q^Q^K^^df^in^p) = qQQ , (89) 

l^K^QPK^^^^,{r^^^fi) = qKQ ■ (90) 

This was the choice adopted in numerical tests verifying the stability of the harmonic IBVP with a plane boundary 
[14| . The powers of r in ([57)l - (|90p are based upon the leading asymptotic behavior of the components for the gauge 
wave ([55)1 and the gravitational waves ([55)1 and (|5S|) . These choices lead to boundary data with the asymptotic 
behavior 



QKK 

IQQ 



qKQ 



f'ju) 



Thus the behavior of qQQ and qQQ imply that the resulting reflection coefficients have overall asymptotic dependence 
no weaker than k = 0{1/R^). 

Our second choice, which is partially suggested by the electromagnetic case ((73|) and leads to weaker reflection, 
consists of the modiflcations 

^K^K^K^'d^ir^ja.p) = QKK , (91) 

iK'^QPK^^d^ir^ja.ff) = qKQ , (92) 

^Q"Q^K^d^{r^jo,p) -1 ^ qQQ, (93) 

Q'^Q'^K^d^jc^p - Q'^K^Q'-d^^a.p - qQQ . (94) 



Now (7.. ~ f{u)/R^ for both gravitational quadrupole waves. For the gauge waves, qQQ ^ f'(u)/R^. Using the Regge- 
Wheeler-Zerilli perturbative formulation and the metric reconstruction method described in [35| we have independently 
checked that this leads to reflection coefficients k — 0{1/R^) for the gravitational waves and k = 0{\/R^) for the 

gauge waves in accord with (|^ . After replacing 7^^ = — /i^^ + ^rj^j.^ and identifying 9^ with the connection Vq of 
the background metric g^f^, (IOT|) - ([M)l correspond to the boundary conditions (|14p - (fT7)) discussed in Sect. II Cl 

Our third choice of boundary conditions, motivated by the first order version of the ^0 boundary condition (|67p . is 



K^^d,,{T^K"K^-l^p) = qKK , (95) 

if^a^(rQ"Q'^7a/3) = qQQ , (96) 

K^'d^irQ^Q^-i^p) = qQQ , (97) 



l^K^d^ir^K'^Q^'j^p) - \Q''d^K^KP^^p + ]^Q^^Q^Q''d^^^p = qKQ 
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However, for the gravitational quadrupole wave (|85p . this leads to qkq ^ f"{u)/B? and so it results in much stronger 
reflection than the first two choices. Thus, as might have been anticipated by the discussion following ([ST|) . the first 
order version of the ^ boundary condition is not as effective as (|94|) -(|9T |) in the case of a spherical boundary. 

V. CONCLUSION 

We have considered the 1 BVP for a coupled system of quasilinear wave equations and established (local in time) well 
posedness for a large class of boundary conditions. In particular, this allows for the formulation of a well posed IBVP 
for quasilinear wave systems in the presence of constraints on finite domains with artificial, nonreflecting boundaries. 
Therefore, we anticipate that our results will have application to a wide range of problems in computational physics. 
Furthermore, since our proof is based on a reduction to a symmetric hyperbolic system with maximal dissipative 
boundary conditions, it also lays the path for constructing stable finite difference discretizations for such systems. 

Our work has been motivated by the importance of the computation of grav itational waves from the inspiral and 
merger of binary black holes, which has enjoyed some recent success [3a, |33, H, |3^, |40] . At present, however, none of 
the simulations of the binary black hole problem have been based upon a well posed IBVP. The closest example is the 
harmonic approach of the Caltech-Cornell group [4l|, [42, Ssf which incorporates the freezing ^q boundary condition 
in second order form and has been shown to be well posed in the generalized sense in the high frequency limit jl5|. 

Our results have potential application to improving the binary black hole simulations. However, many of these 
simulations are carried out using the BSSN formulation [4J, |45J of Einstein's equations, which differs appreciably 
from the harmonic formulation considered here. Although our results constitute a complete analytic treatment of the 
IBVP for the harmonic formulation of Einstein's equations, the extension to the BSSN formulation is not immediately 
evident. For this purpose, it would be useful to reformulate the boundary data for the harmonic problem in terms of 
the intrinsic geometry and extrinsic curvature of the boundary, as has been done for the initial data for the Cauchy 
problem. Such a geometric reformulation remains an outstanding problem. 
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APPENDIX A: SOME USEFUL FORMULAE 

Here we give a short summary of the formulae and conventions underlying the calculational results of See's. IIVBI 
and llVCl We have 



so that 



and 



9a/(«) = -/'(w)if„, u^t-r, K''daKf, = 0. (Al) 



dadp—- = ——K^Kfi + —^{K^rp + r^Kp) + —^r^rp - (— — + -^V^p (A2) 

^.^. . . /(«) f"^^)r, J. 2/>) 6/H 2f(u) 3/(w) 

where Tq, := dar and Vap := dadpr. The spatial components are 

ri — — = (sm 6 cos (f), sin 9 sin <j),cos9), rij — -^ ^ . (A4) 

Our conventions for the polarization dyad give rise to the Cartesian components 

{Q^,Q^, Q^) — (cos6'cos0 — isincj), cos^sin^ + zcos0, — sin0), (A5) 
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which satisfy 

(Q''? + (Qy? ^-sin'^O, Q'^^-Qy-^-ismO, g^- + g^- = sin^f 2cos6'cos<^sin(/) + i(cos2 .^-sin^ </)) ) (A6) 

r r r r \ J 

and 

Q% = ^, Q^a,0^ = ^Q\ Q^d^Q^ = -S^Q^-^. (A7) 

r r r r 

From these follow the necessary commutation relations such as 

[rQf'd^,K''d,]=0. (A8) 
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